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Abstract. We have studied properties and shell structure of the superheavy elements 
from Z — 102 Xo Z — 120 within the framework of the RMF theory. The region 
of study spans nuclides with neutron numbers N — 150 — 190. The Lagrangian 
model NL-SVl with the inclusion of the vector self-coupling of the w-meson has been 
employed in this work. We have performed RMF + BCS calculations for an axially 
deformed configuration of nuclei. The ground-state binding energies, single-particle 
properties and quadrupole deformation of nuclei have been obtained from the mean- 
field minimizations. Two-neutron separation energies, Qa values and a-decay half- 
life have been evaluated. It is shown that a large number of nuclides exhibit the 
phenomenon of shape-coexistence over a significant region of the superheavy elements. 
Shape coexistence of a prolate and an oblate shape is prevalent in nuclides far below 
N = 184, whilst nuclei in the vicinity of = 184 tend to show a shape coexistence 
between a spherical and an oblate shape. The shell structure and 2-neutron separation 
energies obtained with the RMF theory reinforce the neutron number N — 184 as a 
major magic number. It is shown that the neutron number N = 172 acts akin to a 
magic number in the deformed region. It is suggested that the combination Z = 120 
and N = 172 has the potential of being a doubly magic number in the superheavy 
region. 
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1. Introduction 

Significant progress has been made in the discovery of new superheavy nuclei in the 
last decade [Il[2l[3llll[5l[6l[7l[8],|9l |10]. Superheavy elements at the extreme end of 
the periodic table have been synthesized in the laboratory. The synthesis of superheavy 
nuclei with Z = 110-112 at GSI, Darmstadt Pll] and JINR, Dubna P [3 El El [10] has 
stimulated the interest in the properties of superheavy nuclei. Discovery of isotopes of 
Z = 112 and Z = 114 was reported in JINR, Dubna and synthesis of Z = 115 at Berkeley 
was reported some time ago [10]. Earlier superheavy nuclei could be unambiguously 
identified by their a-decay chains leading to already known nuclei. However, the decay 
chains of the newly found superheavy nuclei cannot be linked to any nuclei known due 
to their relatively extreme positions in the periodic table. Therefore, their identification 
depends much upon a comparison with theoretical models. 

Much progress has also been made in theoretical investigations of the properties of 
superheavy nuclei. Theoretical models such as macroscopic-microscopic (MM) models 
have been employed to discern and predict properties of superheavy nuclei. At the 
same time, various mean-field approaches have also been employed to investigate basic 
properties in this domain. The prediction of ?84ll4 was confirmed by MM models such 
as FRDM with Folded Yukawa single-particle potential [11]. For the next doubly- magic 
nucleus, the available self-consistent models predict this to the more proton rich f 72 120 
or the even heavier one at 11^126. Predictions in various models tend to diverge primarily 
due differences in the form of the single-particle potentials. 

Possibly different extrapolations arise in self-consistent Relativistic Mean-Field 
(RMF) theory. Here, the spin-orbit coupling is naturally built-in as a consequence 
of the Dirac-Lorentz structure of nucleons, whereas in all non-RMF models, it has to 
be introduced by hand. The non-RMF models tend to overestimate spin-orbit splitting 
with increase in mass number, thus affecting actual predictions of shell closures in the 
superheavy region. Superheavy nuclei with inherently large density of single-particle 
states can be a sensitive probe for models of nuclear structure. 

Various microscopic approaches such as non-relativistic density-dependent Skyrme 
Hartree-Fock (SHF) theory [El [131 [H m [HI m [18] and that of MM type 
[T9l [20I [211 [22I [23] are used extensively to investigate the properties and structure 
of superheavy nuclei. In the MM approach, a sum of smooth energy based upon liquid- 
drop type formula on which shell correction is imposed using the method of Strutinsky 
is obtained as the total binding energy of nuclei. The Finite-Range Droplet Model 
(FRDM) has achieved a significant success amongst the MM models. The mass formula 
FRDM is used to calculate the shell correction energies in order to identify the major 
magic numbers in the region of superheavy nuclei. A major proton magic number 
Z = 114 and a neutron magic number = 184 have been predicted by the FRDM. 

The MM approach has also been performed by Sobiczewski group [T9 | [20 |[2n[22l[23] . 
where the Strutinsky shell correction [24l [25] is used for the microscopic component, 
whereas Yukawa + exponential model is used for the macroscopic part. Superheavy 
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nuclei have also been studied extensively within the self-consistent mean-field models 
of Skyrme Hartree-Fock (SHF) type. Comparative studies of SHE and RMF models 
have been made. One of the first studies |26] performed within RMF theory earlier on 
was using the force NL-SH [27]. It was shown that neutron number = 184 appears 
to be magic for lighter superheavy nuclei whereas its importance is reduced for heavier 
superheavies with Z > 114. There was also an indication of a small deformed shell- 
closure at Z = 114: for the proton number. In comparison, the FRDM predicted a 
strong shell gap at Z = 114 in the deformed region. 

Extensive shell correction calculations were performed by Kruppa et al. [28] for 
sets of the Skyrme and RMF forces in order to identify the magic numbers in the region 
of superheavy elements theoretically. The Skyrme models predict the strongest shell 
effects for = 184 and Z = 124, 126 and not at Z = 114, whereas RMF theory did not 
show a clear shell gap at = 184 especially for heavier SHE. This was attributed to 
the difference in the shell structure in the above two approaches. This may have been 
due to a different isospin dependence in the Skyrme model as compared to the RMF 
theory [291 EQ] ■ In comparison, all the Skyrme forces have predicted a broad valley of 
shell stabilization around Z = 120 and = 172 and 184. On the other hand, Cwiok et 
al. [31] have suggested a shift in the island of stability towards 184I26. 

In spite of an impressive agreement with experimental data for the heaviest 
elements, the theoretical uncertainities are large when extrapolating to unknown regions 
of the nuclear chart. In particular, there is no consensus among theories with regard 
to the regions of the shell stability for superheavy nuclei. Since in these nuclei the 
single-particle level density is relatively large, small shifts in the position of single- 
particle levels (e.g., due to the Coulomb or spin-orbit interaction) can be crucial for 
determining the shell stability of a nucleus. While most macroscopic-microscopic (non- 
self-consistent) models predict Z = 114 to be magic, most self-consistent calculations 
suggest that the center of the proton shell stability should be moved up to higher proton 
numbers, i.e., Z = 120,124, or 126. For neutrons most non-relativistic calculations 
predict magic gap at A^ = 184 while the relativistic mean-field theory yields A^ = 172 
- due to a slightly different spin-orbit interaction ^6] . The experimental determination 
of superheavy nuclei from this region will thus be very important for pinning down the 
fundamental question of the spin-orbit force. 

A number of theoretical studies have been carried out within the framework of the 
RMF theory [371 EHl [39] to investigate properties of superheavy nuclei, in particular, a- 
decay properties. In ref. [38] the nature of possible magic numbers within the framework 
of spherical relativistic Hartree-Bogoliubov theory using various interactions has been 
studied. A recent theoretical work [17| sheds some light on the question of magic 
superheavy nuclei. According to the calculations, the patterns of single-particle levels 
are significantly modified in superheavy elements. Firstly, the overall level density grows 
with mass number A, as a A^^^. Secondly, no pronounced and uniquely preferred 
energy gaps appear in the spectrum. This shows that shell closures which are to be 
associated with large gaps in the spectrum are not robust in superheavy nuclei. The 
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theory predicts that beyond Z = 82 and = 126 the usual locahzation of shell effects 
at magic numbers is diminished. Instead, theory predicts fairly wide areas of large shell 
stabihzation without magic gaps. 

The Relativistic Mean-Field (RMF) theory [40j has proved to be successful as a 
framework for description of various facets of nuclear properties. In the RMF theory, 
the nuclear force is produced by a virtual exchange of various mesons. The nuclear 
saturation is achieved by a balance between an attractive a- and a repulsive cu-field. The 
relativistic Lorentz covariance of the theory allows an intrinsic spin-orbit interaction 
based upon exchange of a- and cj-mesons. This has been shown to be advantageous 
for properties which depend upon spin-orbit potential. An immediate advantage of 
the proper spin-orbit potential has been the success of the RMF theory to be able to 
describe the anomalous kink in the isotope shifts of Pb nuclei [41]. The isotopic shifts 
of Pb have been measured with a high precision using atomic beam laser spectroscopy 
and are known to show a pronounced kink about the magic number = 126. The 
non-relativistic approaches based upon the Skyrme and Gogny forces have, on the other 
hand, been unable to reproduce this kink. It was shown that this difference in the 
predictions of the Skyrme model and the RMF theory is due to the isospin dependence 
of the spin-orbit term [291 EQ] . The isospin dependence in the RMF theory is provided 
mainly by the coupling constant Qp of the p-meson. However, as the strength of the spin- 
orbit term derives from a large sum of the absolute values of the scalar and vector fields, 
the contribution of p-field is much weaker in comparison. Consequently, the isospin 
dependence of the spin-orbit potential in the RMF theory is relatively weak [29j. With 
the success of the RMF theory in properties of nuclei related to the shell effects, we have 
carried out a study of the superheavy nuclei with an RMF Lagrangian model which has 
been shown to have improved shell effects in nuclei. 

In the present work, we have explored the region of superheavy nuclei within the 
framework of the relativistic mean-field theory using the Lagrangian model with the 
vector self-coupling of the u meson. The region encompasses the superheavy elements 
from Z = 102 up to Z = 120. Nuclides in these isotopic chains with neutron number 
spanning A^ = 150 to A^ = 190 have been considered. We have used the force NL-SVl 
due to Sharma et al. [34j with the vector self-coupling of the u meson. Our focus is 
primarily to look for possible signatures of neutron shell gaps or magic numbers in the 
region of superheavy nuclei. 

2. The Relativistic Mean-Field Theory: Formalism 

In the present work, we have employed the RMF theory to study properties of 
superheavy nuclei. The starting point of the RMF theory is the model due to Walecka 
[lO] , which describes the nucleons as Dirac spinors interacting by the exchange of several 
mesons: scalar mesons [a) couple to the nucleons ('0) through a Yukawa term ipatp and 
produce strong attraction, isoscalar vector mesons (a;) couple to the conserved nucleon 
current ■07/i'0 and cause almost as strong a repulsion. In addition, there are isovector 
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p mesons which couple to the isovector current and photons to produce the well-known 
electromagnetic interaction. 

We restrict ourselves to the Hartree treatment, where the A independent nucleons 
with the single-particle spinors (i= Ir-A) form a Slater determinant and move 
independently in the meson fields. The RMF theory starts with the standard Lagrangian 
density 



+ ^d^ad^a - U{a) - g^tpiO^Ji 

- ^VL^'^VL^^ + ]^mluj^u^ - g^iJiY^^ii^i 

- -^R'^^'R^,^ + ^mjp^p; - gptlJa'^p^f^i 

+ i<?4K^^)^ (1) 

where the summation convention is used and the sum over i runs over all nucleons. 
Isovector quantities are indicated by the arrow bars. As proposed by Boguta and Bodmer 
[12] , the a meson moves in a nonlinear potential of the form: 

U{a) = \n,a'' + ^(72^^ + \gz(y\ (2) 

Here, M, and rup are the nucleon, the a-, the a;-, and the p- meson masses, 

respectively, g^^ guj,gp,a.nd e^/47r=l/137 are the coupling constants for the a-,the u- 
,the p-mesons and for the photon. The coupling constant (74 represents the vector 
self-coupling of a;-meson as introduced in ref. |34j . 

The classical variation principle provides the Klein-Gordon equations for the 
mesons: 

{-A + m2}a(r) = - g^p,{r) - g2a''{r) - gsa^'ir) 

{-A + m'ju^iv) = g^f{v) + g^uj^{v) (3) 
{-A + mJ}p-^(r) =gp3^{v) 
-AA''(r) =ejp^(r), 

whereby the sources are determined by the corresponding density and current 
distributions in the static nucleus. Variation for the nucleonic field gives the Dirac 
equation for nucleons: 

{-la.V + (3M\r) + y(r)}^i(r) = Q^,(r). (4) 

The effective mass M*(r) is determined by the scalar field cr(r) as 

M*{r) = M + g,a{r), (5) 
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and the potential vector potential ^(r) is obtained as: 

V{r) = g^uoir) + g^npoir) + e ~ ""'^ Ao(r). (6) 

Bold faced letters indicate vectors in the 3-dimensional space. These fields are 
the solutions of the inhomogeneous Klein-Gordon equations. As is the case in the 
majority of applications, the contributions of antiparticles are neglected, i.e., the no-sea 
approximation is employed. Equation ([3]) together with Eq. (jl]) provides a closed set 
of equations. In order to describe the ground-state properties of nuclei, static solutions 
are obtained from the equations of motion Eqs. ([3lll]). In this case, the nucleon spinors 
are the eigenvectors of the stationary Dirac equation, which yields the single-particle 
energies e, as eigenvalues. The solution of these equations is obtained iteratively. Using 
these solutions, one calculates physical quantities such as the total energy, charge radii 
and quadrupole moments amomg others. 



2.1. Axially deformed RMF 

For nuclei in the superheavy region exhibit significant deformations, we have employed 
the RMF theory for axially deformed configuration of nuclei. We solve the Dirac 
equation ([3]) as well as the Klein-Gordon equations (jlj) by expansion of the wavefunctions 
into a complete set of eigen solutions of an harmonic oscillator potential |13]. In the 
axially symmetric case the spinors ff" and gf are expanded in terms of the eigenfunctions 
of a deformed axially symmetric oscillator potential 

Vosciz, r^) = ^Muy + ^Mulrl (7) 

Taking volume conservation into account, the two oscillator frequencies hu± and huz 
can be expressed in terms of a deformation parameter (3o. 



The corresponding oscillator length parameters are given by 



and b^ = J— — . (10) 



V V Muj_L 

The volume conservation gives b^^bz = 6g. The basis is now determined by the two 
constants Huq and Po, which are chosen optimally. 

The deformation parameter of the oscillator basis jSo is chosen to be identical for 
the Dirac spinors and the meson fields. The deformation parameter is obtained from 
the calculated quadrupole moments for protons and neutrons through 



Q = Qn + Qp = J^^ARlp2 (11) 
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with Rq = 1.2^4^/'^ (fm). The quadrupole moments are calculated by using the 
expression: 

Qn,p = {2r^P2{cos9))n,p = 2{2z^ -x"- y^)n,p. (12) 

2.2. Shell effect in nuclei 

Shell effects play an important role in identifying and consequently synthesizing 
superheavy elements due to additional stability provided by the shell effects. Hence 
few comments on shell effects would be in order. The shell effects manifest strongly in 
nuclei with the existence of magic numbers. The origin of shell closure has long been 
understood due to spin-orbit coupling and to an ensuing splitting of levels. The spin- 
orbit potential and a bunching of levels create shell closures (magic numbers) which are 
predicted correctly by most of the models. In the non-relativistic density-dependent 
theory of the Skyrme type, the spin-orbit interaction is added phenomenologically and 
its strength is adjusted to reproduce the spin-orbit splitting in ^^O. In the RMF theory, 
on the other hand, the spin-orbit interaction arises naturally due to exchange of a and 
llj mesons between the nucleons. The strength of the interaction is determined by the 
spin-orbit splitting in ^^O and other nuclei, which is uniquely decided by the effective 
mass. As mentioned in the introduction, the shell effects in the RMF theory have been 
advantageous in a few aspects such as isotope shifts in nuclei. In view of this, application 
of the RMF theory to superheavy nuclei may have more realistic predictions. 

2.3. Lagrangian with the vector self-coupling of u-meson 

In this work, we have employed the Lagrangian model with the quartic coupling of 
u meson. The non-linear vector self-coupling of w-meson was introduced by Bodmer 
[B] and properties of nuclear matter were discussed on adding a quartic term in uj- 
meson potential. The inclusion of the vector self-coupling of u meson in addition to the 
non-linear scalar self-coupling of a meson was found to have the effect of softening the 
high-density equation of state (EOS) of the nuclear matter [15]. Subsequently, the force 
NL-SVl [34j with the inclusion of the vector self-coupling of uj meson was developed 
with a view to improve the predictions of the ground-state properties of nuclei, such as 
binding energies, charge radii and isotopes shifts of nuclei along the stability line. It 
was shown that with the Lagrangian model NL-SVl, the shell gaps in nuclei along the 
stability line such as Ni and Sn isotopic chains were improved as compared to Lagrangian 
model with the nonlinear scalar coupling alone [34]. With the improved shell structure 
along the line of /^-stability with NL-SVl, it is expected that it shall have improved 
predictions in the region far away from the stability line. In view of this, we have used 
the parameter set NL-SVl with the vector self-coupling of u meson in our investigations. 
The parameters of the parameter set NL-SVl |34j are given in Table I. 
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Table 1. The Lagrangian parameters of the force NL-SVl [34] used in the RMF 
calculations. 



Parameters 


NL-SVl 


M 


939.0 




510.0349 




783.0 


nip 


763.0 


9a 


10.1248 




12.7266 


9p 


4.4920 


92 


-9.2406 


93 


-15.388 


94 


41.0102 



3. Details of calculations 

The input parameters required to carry out explicit numerical calculations are: neutron 
pairing gap A„, the proton pairing gap Ap and the number of oscillator shells Np and 
of the fermionic wavefunctions and meson fields, respectively. Both the fermionic 
and bosonic wavefunctions have been expanded in a basis of 20 harmonic oscillator shells 
in this work. 

We have carried out a study of superheavy nuclei from No {Z = 102) to Z = 120. 
The neutron number ranges from N = 150 to = 190. This region encompasses the 
neutron number = 184 which has been predicted to be a magic number in various 
theories. A possible deformed shell has also been predicted in this region as discussed in 
Section 1. We have carried out deformed RMF+BCS calculations for superheavy nuclei 
in the above region. For pairing gaps, we have used the formula due to MoUer and Nix 
[i6] as given by: 

An = 4.8 Ar-i/3 (^^3^ 
Ap = 4.8 Z-i/^ (14) 

An axially symmetric deformed configuration with refiection symmetry has been 
assumed for nuclei. For each nucleus an initial basis deformation of prolate as well 
as oblate type has been taken in order to seek a minimum both in the prolate as well as 
in the oblate region of the deformation space. Consequently, we have obtained a prolate 
and an oblate minimum for most of the nuclei. However, nuclei close to a potential 
magic number assume a spherical shape. 
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4. Results and Discussions 

4.I. Binding energies and deformations 

The RMF+BCS minimizations with NL-SVl give rise to a prolate and an oblate solution 
for most of the nuclei. We show the quadrupole deformation p2 for the lowest energy 
state (ground state) for all the isotopic chains in Figures [D^a)- [T](j). A large number of 
nuclei in these chains assume a prolate shape in the ground state. This is especially the 
case in the region of neutron number ~ 150 — 170. Typically, these nuclei exhibit a 
prolate deformation of (32 ~ 0.2 — 0.3. 

For the isotopic chains of No {Z = 102), Rf (Z = 104) and Sg {Z = 106), 
there is a transition from a prolate shape to an oblate shape at ~ 170. Beyond 
= 172, nuclei assume an oblate shape in the ground state for these isotopic chains 
(see Figs. [11(a)— [H^d)). The magnitude of the oblate deformation decreases as the 
neutron number increases above = 172, which approaches a vanishing value as the 
neutron number A^ = 184 arrives. In the vicinity of A^ = 184, nuclei exhibit a spherical 
shape. Such a feature has been commonly predicted by several relativistic as well as 
non-relativistic theories, whereby A^ = 184 constitutes a major magic number in the 
region of superheavy nuclei. 

As one goes to the heavier elements above Sg {Z = 106), one can visualize that the 
tendency of prolate shapes in the range of neutron numbers A^ = 150 — 176 continues 
for heavier superheavy elements as well. The point of shape transition from prolate to 
oblate seems to be shifting to a higher neutron number. For instance, for Hs [Z = 108) 
the prolate to oblate transition occurs at A^ = 178. For the isotopic chain Ds {Z = 110) 
the shape transition point advances to A^ = 180. However, for Z = 112 and above, the 
transition point falls back to A^ = 178 and below it. For instance, for Z = 116 and 
Z = 118, the shape transition occurs at A^ = 174 and A^ = 172, respectively. Similarly, 
the shape transition point for Z = 120 falls back to A^ = 170. Thus, the number of nuclei 
which exhibit an oblate shape below the neutron number A^ = 184 increases steadily in 
going from Z = 112 to Z = 120. As far as the spherical shape for A^ = 184 nuclei is 
concerned, it is exhibited by all the isotopic chains from Z = 102 to Z = 120. This is 
reminiscent of a magic shell character for the neutron number 184. 

The general feature of shape transition from prolate to oblate in the region of 
A^ = 172 — 176 prevails for all the isotopic chains considered in this study. We have 
made a comparison of NL-SVl deformations with those predicted by the macroscopic- 
microscopic mass formula FRDM. For the isotopic chains Z = 102 — 112, the FRDM 
exhibits a predominantly prolate shape for the neutron numbers A^ = 150 — 176. This 
feature is very similar to that predicted by NL-SVl. 

The FRDM shows a shape transition from prolate to an oblate shape at the neutron 
number A^ ~ 178 for the chains Z = 102 — 110. There is only a single nuclide exhibiting 
an oblate shape in FRDM at A^ ~ 178 for these chains, beyond which all the nuclei 
assume a spherical shape for the isotopic chains of Z = 102 — 112 in the FRDM 
predictions, in contrast to NL-SVl. The oblate deformation for A^ = 178 nuclei in this 
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Figure 1. The j32 values for isotopes of the elements Z = 102 — 120 obtained with 
NL-SVl. The results with FRDM are shown for comparison. 



region amounts to ~ 0.1 in the FRDM. This is unhke NL-SVl, where a significant 
number of nuclei exhibit a moderately deformed oblate shape until the supposedly magic 
number = 184 arrived at. 

For the isotopic chains Z — 114 — 118, the behaviour of FRDM for deformation 
properties is rather uneven. For the region of the lower neutron numbers, some nuclei 
in these chains show a moderately prolate shape, whereas a few others exhibit a small 
prolate deformation. This is followed by a shape transition to an oblate deformation 



Superheavy nuclei with the vector self- coupling of the u-meson in the relativistic mean-field theoryll 

at N = 172, followed by a spherical shape in going to higher neutron numbers. That 
nuclei in the vicinity of = 184 are predominantly spherical in FRDM points to the 
hypothesis that = 184 maintains its magicity even for heavier superheavy elements 
in FRDM. 

4-2. Comparison with ETF-SI results 

The mass formula Extended Thomas-Fermi with Strutinsky integral (ETF-SI) [17] 
was obtained to fit experimental masses of more than 1600 nuclei from all over 
the periodic table. The ETF-SI, due to its partly microscopic structure, has found 
broader acceptance for nuclear physics applications. It predicts binding energies and 
deformations of nuclei over a broad range of masses quite well. Here we compare 
the ground state deformations obtained with NL-SVl with those predicted by ETF- 
SI. The quadrupole deformations (32 are compared for the given isotopic chains in 
Figs. [2](a)-l2]^j). Figs. [2]^a)-l2](d) show that ETF-SI deformations are in good agreement 
with the NL-SVl values from = 150 to ~ 166 for the elements from Ds {Z = 102) 
to Hs (Z = 108). For Ds {Z = 110), the agreement between NL-SVl and ETF-SI is 
good up to A^ = 164 (see Fig. Wie)). 

A constrasting picture with the ETF-SI is that for neutron numbers A^ ~ 166 
and above, it shows a strong divergence from the NL-SVl predictions. ETF-SI does 
not exhibit a prolate to oblate transition as do NL-SVl and FRDM predictions. In 
constrast, ETF-SI predicts a constant and a large value of ~ 0.4 for nuclei in the 
range of A^ ~ 168 — 180 for the isotopic chains of Z = 102 — 110. For the heavier 
elements Z = 112 — 120, ETF-SI predicts a nearly constant value /?2 ~ 0.4 for all the 
neutron numbers A^ = 150 — 180. 

ETF-SI predicts a spherical shape of nuclei in the vicinity of A^ = 184 for the 
elements Z = 102 — 114. This feature of ETF-SI is similar to that of FRDM. However, 
for Z = 116, ETF-SI bypasses the spherical shape for A^ = 184. On the other hand, for 
several isotopic chains ETF-SI predicts a shape transition from spherical to a strongly 
prolate shape at A^ = 188 - 190 (see Figsl2](b)- EJ^g)). 

For Z = 110 and Z = 112, nuclei with two neutrons below and two neutrons above 
A^ = 184 exhibit a spherical shape, whereas for Z = 114 there is only one nuclide i.e, 
A^ = 184 that shows a spherical shape. It's neighbours on both sides, assume a prolate 
shape in ETF-SI. This is indicative of a decrease in the shell strength at A^ = 184 with 
ETF-SI. This point is supported by a deformed shape taken by the A^ = 184 nuclide for 
Z = 116 (See Fig. El^h).) with ETF-SI. Whilst ETF-SI values for Z = 118 and Z = 120 
are not known, the aforesaid diminishing of A^ = 184 shell strength points to an erosion 
of a possible A^ = 184 magicity for the heaviest superheavy elements from Z = 116 to 
Z = 120 with ETF-SI. 

A summary of quadrupole deformation in the ground state of superheavy nuclei 
with NL-SVl obtained in this work is shown in the FigO The nuclei denoted by the 
red colour indicate a prolate ground state. An oblate ground state is shown by the 
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Figure 2. The P2 values for the isotopes of the elements Z = 102 — 120 obtained with 
NL-SVl are compared with ETF-SI predictions. 



green colour. Spherical nuclei are indicated by yellow. As seen in Fig. [3l nuclides with 
A'^ = 150 — 172 for almost all the isotopic chains exhibit a prolate shape. On the other 
hand, nuclides with > 172 exhibit an oblate shape until a spherical shape appears for 
nuclei close to = 184. There is a little exception here in that for SHE with Z = 108 
to Z = 112, a few nuclides above A = 172 continue to exhibit a prolate shape. This 
indicates an intrusion of prolate shape into the domain of oblate shapes as indicated in 
the middle of the figure. Nuclei in the vicinity of A^ = 184 are predominantly spherical 
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Figure 3. The landscape of the quadrupole deformation P2 obtained with NL-SVl 
for Z = 102- 120 and TV = 150 - 190. 



as shown by the yellow blocks. In a few cases and especially for Z = 102, nuclei two 
neutrons below and two neutrons above = 184 assume a spherical shape. As shown 
by the picture, = 184 acts as a major magic number in the neutron-rich region for 
all the SHE with Z = 102 - 120. 

As we will see in the single-particle levels, the neutron number = 184 remains 
as a major magic number. Also we will see that in the RMF theory with NL-SVl, 
the proton number Z = 120 is also a magic number. The combination Z = 120 and 
= 184 seems to act as double closed shell in the RMF theory. 

Nuclides above A^ = 184 show a prolate (red) region for a few lighter SHE. However, 
for heavy SHE with Z = 108, nuclides above A^ = 184 are predominantly oblate (green) 
as shown in the Fig. [31 Thus, a significant region below as well as above A^ = 184 lends 
itself to an oblate deformation. 



4-3. Shape co-existence in the superheavy region 

The phenomenon of shape co-existence in the region of deformed nuclei arises due to 
interplay of deformed single-particle levels. Splittings and alignment of deformed single- 
particle levels leads often to ground states with a prolate and an oblate shape nearly 
degenerate in energy. For practical purposes, when the two minimum energy states are 
within ~1 MeV, the two states are said to be shape co-existent. 

We show in Figs. Hl^a)— |l](j) nuclides which exhibit a shape co-existence with NL- 
SVl. Looking at the figures one can see that a large number of nuclides exhibit the shape 
co-existence for the superheavy chains. The nuclides in isotopic chains from Z = 102 
(No) to Z = 112 show a common pattern of shape co-existence between prolate and 
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Figure 4. Isotopes of Z = 102—120 exhibiting the shape-coexistence in the Lagrangian 
model NL-SVl. 



oblate shapes for the region of neutron numbers ~ 170 to ~ 182. Whilst for 
Z = 102 and Z = 104 shape co-existence is limited to neutron numbers ~ 168 to 
~ 178, for Z = 106 to Z = 112 the region of shape co-existence extends to neutron 
number A^ = 182. Whereas for the lighter neutron numbers there is a prolate-oblate 
shape co-existence for elements Z = 102 to Z = 112, there is a shape co-existence of 
spherical-oblate shapes for A^ = 182 isotones from Z = 104 — 112. This indicates a 
softer nature of A^ = 184 shell closure even for the lighter superheavy elements. Several 
nuclides with A^ > 184, also exhibit a spherical-oblate shape co-existence for the isotopic 
chains Z = 106 - 112. 

Figures 111(g)- ll](j) show the shape co-existence for heavier superheavy elements from 
Z = 114 to Z = 120. Unlike the picture for superheavy elements Z = 102 — 112, the 
shape co-existence in the lighter neutron region is reduced significantly for Z = 114—120. 
Only a few nuclides below A^ = 180 exhibit a prolate-oblate shape co-existence. For 
Z = 114 and Z = 116, several nuclides exhibit a spherical-oblate shape co-existence 
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Figure 5. The landscape of shape co-existence in superheavy nuclei with NL-SVl. 
The prolate-oblate and spherical-oblate shape co-existenec is depicted. 



for neutron numbers > 184. Interestingly, = 184 nuclide for Z = 116 shows a 
spherical-oblate shape co-existence. This indicates a softening of the shell effects for 
= 184 for Z = 116. 

The superheavy elements Z = 118 and Z = 120 exhibit the shape co-existence 
rather scarcely i.e, only two isotopes = 170 and A^ = 172 show a prolate-oblate shape 
co-existence. In the rest of the neutron region, there is no shape co-existence for the 
heaviest superheavy elements with Z = 118 and Z = 120 This may imply that the shell 
strength at A^ = 184 has picked up for these chains and that a readjustment of the 
deformed single-particle levels does not lend to degenerate shapes in the ground state. 

A summary of the shape co-existence in superheavy nuclei is illustrated in Fig. [5l 
The brown squares depict a shape co-existence between a prolate and an oblate shape, 
whereas the blue squares represent a shape co-existence between a spherical and an 
oblate shape. As the figure shows, the phenomenon of prolate-oblate shape co-existence 
pervades the region from A^ = 170 to A^ = 180. For the lighter elements this shape co- 
existence lies in the neutron region 170-178. The lower boundary of this region gradually 
shifts to higher neutron numbers Z = 110 is reached. As one moves above Z = 110, the 
region shifts gradually towards lower neutron numbers again. From Z = 114 to Z = 120, 
however, there are fewer cases of the prolate-oblate shape co-existence. Apparently, the 
region of prolate-oblate shape co-existence as shown by the brown colour is symmetric 
about Z = 110. In addition, there are a few more cases of prolate-oblate shape co- 
existence for higher neutron numbers near N= 190. 

The region of spherical-oblate shape co-existence lies near A^ = 184. The neutron 
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number = 184 being a major magic number, the ground-state of nuclides with 
= 184 is spherical as seen earlier. For almost all the isotopic chains there is no 
shape co-existence (empty squares) at = 184, with the exception of Z = 116. Only 
for Z = 116, there is a spherical-oblate shape co-existence at = 184. This is due to 
an erosion of the magicity of A^ = 184 for Z = 116. 

Nuclides with A^ = 182 show a spherical-oblate shape co-existence as shown by 
blue squares. Nuclides with A^ = 186 and A^ = 188 also exhibit a spherical-oblate shape 
co-existence for several isotopic chains. The fact that nuclides just below and above 
A^ = 184 display a spherical shape co-existence indicates that A^ = 184 is not such a 
strong magic number. This can be constrasted with other magic numbers in the periodic 
table which exhibit robustness. There are rarely occasions of a shape-coexistence or even 
a deformation in the direct vicinity of a major magic number especially when one is not 
very far from the line of /3-stability. Superheavy nuclei especially those associated to 
A^ = 184 do not fall in this category and are deemed to be in a region far from the 
stability line. Therefore, mellowing of the shell strength of A^ = 184 is not unexpected. 

4-4- Two-neutron separation energies 

Magic numbers in nuclei are characterized by a large shell gap between the last single- 
particle level of the magic number and the next single-particle level above it. This results 
in a larger value of 2-neutron separation energy for a nucleus with 2-neutrons more 
than the magic number. Consequently, the major magic numbers exhibit a characteristic 
kink in S2n values all over the periodic table. Thus, the difference S2n{N) — S2n{N -\- 2) 
at the magic neutron number A^ reflects its shell gap. Its magnitude is a measure of the 
shell strength. 

With a view to visualize as to whether there are possible shell closures in neutron 
numbers for superheavy nuclei, we have computed S2n values for all the isotopic chains. 
Using the binding energies of neighbouring isotopes, 82™ is defined as, 

S2n = B.E. (N)- B.E {N + 2). (15) 

The S2n values for the isotopic chains considered in this work are shown in Fig. El The 
values correspond to the lowest energy state (ground state) obtained with RMF-I-BCS 
minimization using the force NL-SVl. As expected, the S2n values show a decreasing 
trend in going to nuclei with higher neutrons. The S2n values show a rather monotonous 
decrease in going from A^ = 152 to A^ ~ 170 with a slight kink-like structure in the 
neighbourhood of A^ ~ 154 and A^ ~ 164. Due to significant changes in deformation 
properties of nuclei near A^ ~ 172 and A^ ~ 184, there is much more structure in the 
S2n values. 

In order to visualize the structure in the S2n values, we have split Fig. [16] into two 
parts which are shown in Fig. [71 In Fig. [T7| (a) the S2n values for A^ = 152 — 168 are 
shown. The higher neutron part with A^ = 168 — 190 is shown in Fig. [71 (b). One 
can notice a slight kink at A^ = 154 for several isotopic chains. A kink like structure 
is visible for Z = 106 — 116. This structure diminishes for Z = 118 and Z = 120. 
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The structure at = 154 is compatible with a deformed shell gap predicted by non- 
relativistic macroscopic-microscopic approaches [T3] . 

One can also notice a slight kink at ~ 164. Thus, there seems to exist a few 
pockets of presumably deformed shell gaps at ~ 154 and ~ 164. However, there 
are no signatures of any strong magic number in the region A^ = 152 — 168. Fig. [7] 
(b) shows structures in S2n values at A^ ~ 172, 178 and 184. The behaviour of S2n 
values in the neutron region A^ = 168 — 190 is far from monotonous. However, for the 
lighter superheavy elements, there is very little structure at A^ ~ 172 except for Z = 102 
and 104. For the heavier elements with Z = 116 — 120, the kink at A^ ~ 172 shows a 
slightly upward trend. A kink like structure in this region indicates a possible A^ ~ 172 
deformed shell gap for superheavy elements. A combination of A^ ~ 172 and Z ~ 120 
might lend itself to a structure akin to a double shell closure. 

The region near A^ ~ 184 is full of structural effects as can be seen in Fig. [7| (b). 
A kink in the S2n values can be seen clearly at A^ = 184 in several isotopic chains with 
Z = 104 — 110. Thus, for the lighter superheavy elements NL-SVl predicts a shell 
closure at A^ = 184 as indicated by the kink. For the deformed RMF+BCS calculations 
with NL-SVl, a sharp decrease in the values at A^ = 184 is indicative of a shell 
closure which is consistent with that with NL-SH [26] . For SHE above Z >112 the kink 
and the structure at A^ = 184 is reduced significantly. Thus, for the heaviest SHE, 
values do not give an unambiguous indication of a shell closure at A^ = 184 . 
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Figure 7. (a) Two- neutron separation energy S2n for the isotopes of elements Z = 102 
to Z = 120, with neutron numbers ranging from N = 152 to = 170. (b) 82^ values 
for = 170 to = 190. 



4-5. Alpha- decay properties - Qa values 

In the pursuit of superheavy heavy elements, it is one of aims of research to find out 
regions of extra-stabihty with a view to synthesize superheavy elements with the most 
appropriate projectile and target combinations. With this in mind, it will be helpful if 
one can explore the regions of higher stability theoretically. Depending upon the region 
of an extra-stability which would arise from shell gaps and possible shell closures, the 
half-life of a-decay is a potential indicator of a possible value of stability for an area 
of enhanced stability. It is expected that for nuclei near a shell closure the ct-decay 
half-lives would be larger than their neighbours. 

Nuclei in the transuranium and superheavy region are known to be a emitters. 
The rate of a-decay depends strongly on Qa value of a nucleus. We have calculated 
Qa values for nuclei in all the isotopic chains we have considered. The values thus 
obtained are shown in Fig. [HI As one moves from Z = 104 to Z = 120, Qa values show a 
continuous increase with increase in the charge number of a nucleus. For all the isotopic 
chains one can see a decrease in the Qa values with an increase in neutron number. A 
decreasing a-value implies a higher ct-decay half-life. Consequently, in moving from the 
lower neutron number = 150 towards = 184, one has a stabilizing factor with a 
decrease in the Qa values. For nuchdes with A^ = 186, a slight peak in Q^-values for 
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Figure 8. Qa values for the isotopic chains with Z = 102 — 120 and N = 150 — 190 
obtained with NL-SVl. 

several isotopic chains with Z — 104 — 114 can be seen. In terms of a-decay half-hfe 
(as wc will see below) it will mean that nuclides with N — 186 become more unstable 
towards a-decay. 

In going above = 184 one sees a little peak at = 186 for the isotopic chains 
with Z = 106 — 114, though for the chain with Z = 104, Qa shows an increasing trend 
from = 184 to = 186. In terms of a-decay half-life this would imply a decreased 
half-hfe for nuchdes with N — 186. For the isotopic chains above Z — llA such a 
behaviour is not seen, but for only Z — 120. Above N — 186, there is a decrease in Qa 
values for all the isotopic chains. 

4-6. Comparison with the experimental data (systematics) 

A large amount of experimental effort to the synthesis of superheavy elements as well 
as to explore the nuclear properties of superheavy nuclei by the Lawrence Berkeley 
National Laboratory, Berkeley, GSI, Darmstadt and the Joint Institute for Nuclear 
Research, Dubna. Binding energies (masses), Qa values and a-dccay half-lives of many 
nuclei have been determined or estimated. Naturally, such data is accessible only for 
light mass SHE especially with Z — 100 — 110. Only few data exist for Z = 112 and 
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Figure 9. values obtained with NL-SVl for the isotopic chains with Z = 104, 106 
and 108 for the neutron numbers ranging from N = 150 — 190. The experimental 
values taken from the systematics of Audi-Wapstra [15] are shown for comparison. 



We show in Figs. [H] and (TUl Qa values obtained with NL-SVl. These are compared 
with Qa values taken from 2003 compilation of systematics by Audi-Wapstra [H]. As 
seen in Fig. [HI NL-SVl values underestimate the empirical Qa values for Z = 104 — 108 
slightly. The disagreement was a bit more pronounced for Z = 104. For Z = 106 the 
agreement of NL-SVl values with the data is reasonably good. For the case oi Z = 104 
and Z = 106, the kink in Qa values at = 162 in the experimental data is reproduced 
well by NL-SVl. For the chain of Z = 108, though NL-SVl values are close to the 
experimental data, the kink at A = 162 in the experimental data is not reproduced 
by NL-SVl. This may be due to complex shell structure around Z = 108 than that 
presented by the mean field in the RMF. 

In Fig. [To] the NL-SVl values show agreement with several data points, whereas 
there is a discrepancy of more than 1 MeV for the data points of A = 158 and A = 166. 
Such differences may arise from a different evolution of shape in reality than predicted 
by the theory. In contrast, the NL-SVl values show a very good agreement with 
experimental data for SHE chain with Z = 112. It may however be mentioned that 
for Z = 112 the empirical data available is limited to four data points. The number of 
data points available for Z = 114 is limited to two points only. The trend of these two 
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data points is represented quite well by the theoretical values, though the discrepancy 
for data point at = 174 is little more than 1 MeV. Given the complex picture of 
superheavy region that evolves due to interplay between deformation and single-particle 
states, the theoretical description of values as seen in Figs.[9]and [TO]can be considered 
as satisfactory. 
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Figure 10. Qa values obtained witii NL-SVl for the isotopic ciiains oi Z = 110, 112 
and 114 for tlie neutron numbers N = 150 — 190. Tlie experimental values taken from 
the systematics compiled by Audi-Wapstra [48] are shown for comparison. 



4-7. Alpha- decay half-lives 

Using the Qo-values, a-decay half-lives can be calculated using the phenomenological 
formula of Viola and Seaborg as given by, 

logT^ = {aZ + b)Q-^^'^ + {cZ + d), (16) 

where Z is the atomic number of the parent nucleus and Qa is the a-decay energy in 
MeV. Then, is obtained in seconds. The parameters a, b, c and d have been fitted to 
a-decay half-lives of known nuclei over a large range. This includes new experimental 
data obtained in the last decade. The values of a, b, c and d have been obtained by 
Sobiczewski et al. [IH]. The values of these parameter are a = 1.66175, b = —8.5166, 
c = -0.20228, and d = -33.9069. 
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Using the Viola and Seaborg formula of Eq. (fTB]l , we have calculated the a-decay 
half-lives T„ for nuclides in all the isotopic chains. The results are shown in Fig. [TTJ 
As can be seen from Eq. ( |T6l) . log Ta is inversely proportional to \JQa- The decreasing 
trend of Qa with neutron number as shown in Fig. [8] translates into an increasing trend 
for log Ta. As seen in Fig. [TTl log values show a rather linear increase with neutron 
number in the region = 150 to A^ = 180 for all the isotopic chains with the exception 
of the heaviest SHE Z = 116 — 120. The latter show a linear behaviour in the curves in 
the neutron region A^ = 150 — 170. Above = 170 the slope exhibits a further increase 
in the value of log T^. A nearly linear increase of log Ta with neutron number implies 
an exponential increase in Ta in going to higher neutron numbers. The half-life curve 
shows the highest values for the isotopic chain Z = 102. As the Z value of an isotopic 
chain increases, the half-life curves show a gradual decrease. This shows that in going 
to heavier SHE, nuclides become less and less unstable. 
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Figure 11. The log Ta values for the isotopic chains with Z = 102 — 120 and 
iV = 150 - 190 obtained with NL-SVl. 

Commensurate with values in Fig. [8], one notices a few peak-like structures in 
half-lives at certain neutron numbers. A small peak-like structure can be seen at the 
neutron number A^ = 154 in nearly all the isotopic chains with the exception of the 
heaviest chains Z = 118 and Z = 120. The increased Tq, at A^ = 154 can be understood 
in terms of a deformed shell gap at this neutron number. A similar behaviour (feature) 
has also been predicted by non- relativistic macroscopic- microscopic calculations |50j . 
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Another minor enhancement in can be seen in the vicinity of ~ 160 — 162 
for the isotopic chains for the hghtest SHE Z = 102 — 112. The region of ~ 162 has 
also been predicted to be a deformed shell ^J. A major enhancement in Ta values can 
be conspicuously seen near = 182 — 184. A peak like structure in log curves can 
be seen near ~ 184 for the isotopic chains of SHE with Z = 102 — 112. Apparently, 
A^ = 182 exhibits a higher value of as compared to A^ = 184. In going to A^ = 186, 
there is a sudden decrease in values by a few orders of magnitude especially for lighter 
SHE with Z = 104 — 108. For the chains with Z = 110 — 112, there is also a decrease 
in Ta a.t N = 186. This decrease is, however, smaller than for lighter SHE. Thus, for 
SHE with Z = 102 — 114, the region of neutron number A^ = 184 denotes a region of an 
enhanced stabihty against a-decay. 

As one goes above Z = 114 to heavier SHE, the Tq values show a significant decrease 
by several orders of magnitude. Nonetheless, SHE with Z = 116 to Z = 120 still show 
a maximum near A^ ~ 184. This implies that even for the heavier SHE the region of 
A^ ~ 184 shows a enhanced stability as compared to their lighter neutron counterparts. 
Thus, the region with A^ = 184 would be suitable for synthesizing heavier superheavy 
elements. For the case of Z = 120 isotopic chain, there is also a peak in the Ta curve 
in the Fig. [11] at A^ = 174. This indicates that the combination Z = 120 and A^ = 174 
may be suitable for synthesizing the Z = 120 superheavy element. 

4-8. Potential energy landscapes 

A potential energy landscape represents the total energy of a nucleus as a function 
of deformation. It reflects how a nucleus rearranges its single-particle structure under 
the influence of deformation. We have carried out constrained RMF calculations with 
the Lagrangian set NL-SVl. The quadrupole deformation is used as a constraint 
parameter i.e, the binding energy of a nucleus is calculated for a given deformation. 
The RMF calculations are carried out for each value of /?2 in the deformation plane with 
a small grid in (32- 

We show in Fig. [12] the potential energy surface (PES) for a nuclide of the lighter 
superheavy element lo^Rf {Z = 104) with the neutron number A^ = 174. Two valleys 
can be seen, one in an oblate region and another in the prolate region. This picture 
is a clear reflection of a prolate-oblate shape co-existence whereby the two minima are 
nearly degenerate in energy. The actual ground state in this case is, however, prolate 
as shown by the slightly lower valley in the prolate region. The PES for a heavier SHE 
for the nuclide ^^'^116 {Z = 116) with the neutron number A^ = 184 is shown in Fig. [13] 
This figure illustrates a shape co-existence between a spherical and an oblate shape, 
as has been indicated by a blue square in Fig. [5] However, here the ground state is 
spherical in shape. This is consistent with A^ = 184 still retaining some magicity in 
going to Z = 116. 

We show in Fig.[ll]the potential energy landscape for another nuclide with Z = 118 
and A^ = 172. It shows a clear minimum for a spherical shape (/32=0). There is no shape 
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Figure 12. The potential energy landscape for ^^^Rfi74 obtained with NL-SVl. 




> 
CD 



< 
UJ 



-7.08 



-7.09 - 




Figure 13. The potential energy landscape for "^""116184 obtained with NL-SVl 



co-existence with any other shape for this nucleus. The spherical shape aX N — 172 
reflects the character of a shell closure in the midst of deformed nuclei in a very heavy 
superheavy region. A spherical shape for = 172 in the deformed region is accentuated 
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Figure 14. The potential energy landscape for ^^°118i72 obtained with NL-SVl. 

by the proximity of its proton number to the presumed magic number Z = 120. 
4-9. The single-particle levels 

The single-particle structure of nuclei is a good reflection of magic numbers, whereby 
larger shell gaps at a magic number can be seen in single-particle levels. In order to 
visualize some of the potential magic numbers in the region of superheavy nuclei, here we 
show the single-particle spectrum for a few key nuclei. As mentioned before, the neutron 
number = 172 is perceived as akin to a magic number. We show in Fig. [15] neutron 
single-particle levels calculated for nuclides with N = 172 for SHEs with Z = 112 — 120. 
In all these nuclei the shell gaps at = 172 and at A = 184 can be seen clearly. In 
these nuclei the shell gap at A = 172 is maintained in going from Z = 112 to Z = 120. 
Quantitatively, there is no change as such in the shell gap at = 172 in going from 
Z = 112 to Z = 120. However, the gap itself is pushed down in energy in reaching 
Z = 120. Thus, the neutron number A^ = 172 would be helpful in synthesizing the 
superheavy elements with Z = 118 and Z = 120. 

In Fig. [16] we show the neutron single-particle spectrum for nuclides with A = 184 
for SHEs with Z = 112 — 120. The shell gap at A = 184 can be seen clearly. Thus, 
A = 184 manifests as a major magic number. This is consistent with several other 
theoretical predictions. In the RMF theory with NL-SVl, the shell gap at = 184 
shows a decrease in going from Z=112toZ=120. We also see a major shell gap at 
A^ = 198. However, it may be of little practical importance in synthesis of superheavy 
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Figure 15. The neutron single-particle levels for N ~ 172 isotones for Z = 112, 114, 
116, 118, 120 with NL-SVl. 
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Figure 16. The neutron single-particle energy levels for Z = 112, 114, 116, 118, 120 
and N = 184 with NL-SVl. 



elements due to very large number of neutrons involved. 

In order to visualize the single-particle structure of SHE Z = 120, we show 
in Fig. [T71 the proton single-particle levels obtained with NL-SVl for nuclides with 
N = 172, N = 184 and = 198. The nuclide with Z = 120 and N = 172 is very close 
to being unbound. The proton shell gap at Z = 120 can be seen for = 172, though 
it couples to continuum. From a practical point of view, the nuclide with Z = 120 
and = 172, which shows the character of a doubly magic nucleus can be useful in 
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Figure 17. The proton single-particle levels for the isotopes of Z = 120 with N = 
172, 184 and 198 using NL-SVl. 

synthesizing SHE with Z = 118 and Z = 120. The combination Z = 120 and = 184 
also comes out as a doubly magic configuration. It may be possibly of use in synthesizing 
the heaviest superheavy nuclei. 

5. Summary and Conclusions 

We have studied the structure of superheavy nuclei within the framework of the 
relativistic mean-field theory. The region of superheavy elements from Z = 102 
to Z = 120 has been explored. The isotopic chains of the superheavy elements 
encompassing the neutron number from = 150 to = 190 have been investigated. 
The Lagrangian model NL-SVl with the inclusion of the vector self-coupling of the 
a;-meson has been employed. 

RMF calculations for a large number of nuchdes in the isotopic chains of superheavy 
elements have been performed by taking an axially symmetric deformed configuration. 
The binding energies, deformation properties and single-particle levels have been 
obtained from RMF+BCS minimizations. Nuclides in general are found to possess 
a moderate oblate or a prolate deformation. On the other hand, nuclei in the vicinity 
of the neutron number = 184 display a spherical shape. Display of a spherical shape 
for all the isotopic chains indicates an inclination of A^ = 184 being toward a potential 
magicity. 

It is shown that nuclei in the region of superheavy elements exhibit the phenomenon 
of shape co-existence. A large number of nuclei much below N — ISA exhibit a shape 
co-existence between a prolate and an oblate shape. This is specially the case for the 
isotopic chains from Z — 102 to Z — 114. Some isotopic chains exhibit a shape- 
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coexistence between a spherical and an oblate shape especially those in the vicinity of 
N = 184. The shape-coexistence in nuclei is demonstrated by the potential energy 
landscapes obtained for a few nuclei. 

The Qa values and the corresponding a-decay half-lives were obtained using the 
ground-state binding energies. The a-decay half-lives Tq, demonstrate clearly that in 
the vicinity of = 184 nuclides have a significantly larger as compared to their 
neighbours. This is again a strong indicator of the magicity of the neutron number 
N = 184. 

Two-neutron separation energies of nuclides were calculated using the results with 
NL-SVl. The S2n values exhibit a small kink near = 160 and = 172. This 
suggests that the neutron number = 160 and A^ = 172 have a larger shell gap than 
their neighbours. In comparison, a strong kink in S2n values at N — 184 is a clearest 
indication of a major shell gap at A^ = 184. 

The single-particle levels obtained with the Lagrangian model NL-SVl show that 
there exist shell gaps at neutron numbers A^ = 172 and N — 184. Whilst A^ = 172 can 
not be construed as a major magic number as displayed by a mild kink in values, 
the neutron number N — 184 does behave as a major magic number in the neutron- 
rich region. Evidently, the magicity oi N — 184 is demonstrated succintly by a large 
shell gap at A^ = 184 seen in the single-particle levels. This picture of the magicity of 
N=184 in the RMF theory is consistent with the predictions in the density-dependent 
Skyrmc Hartree-Fock approach and also with the previous studies performed within the 
RMF theory using the conventional model of the nonlinear scalar self-couplings. Thus, 
our study with the Lagrangian model NL-SVl with the vector self-coupling of u meson 
reinforces the magic character of the neutron number N — 184. Albeit, a decrease in 
the shell gap a,t N — 184 in going towards Z — 120 imphes a weakening of the shell 
strength in approaching the extreme end of the periodic table. 

Whilst a combination of N — 184 in conjuction with Z in the vicinity of 120 pro- 
duces nuclei which have impractically large number of neutrons, the combination of the 
neutron number A^ = 172 and Z = 120 does lie within the domain of feasibility. Our 
study shows that this combination is suggestive of a double magicity. Thus, synthesis 
of superheavy elements heavier than Z = 116 with neutron number close to A^ = 172 
should be susceptible to an extra stability rendered by the apparent magic nature of 
Z = 120 and A^ = 172. 
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